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Tableaux combinatorics of the two-species PASEP 

Olya Mandelshtam and Xavier Viennot 

Abstract 

We study a two-species PASEP, in which there are two types of particles, “heavy” and “light”, 
hopping right and left on a one-dimensional lattice of n cells with open boundaries. In this 
process, only the “heavy” particles can enter on the left of the lattice and exit from the right 
of the lattice. In the bulk, any transition where a heavier particle type swaps places with an 
adjacent lighter particle type is possible. We generalize combinatorial results of Corteel and 
Williams for the ordinary PASEP by defining a combinatorial object which we call a rhombic 
alternative tableau that gives a combinatorial formula for the stationary probabilities for the 
states of this two-species PASEP. 


1 Introduction 

The asymmetric exclusion process (ASEP) is a model of interacting particles introduced in biology 
and mathematics in the late 1960’s mm- The particles may hop to the right and left in a hnite lattice 
of cells and enter and exit at the boundaries of the lattice, with the exclusion condition which means 
there can be at most one particle in each cell. In the most general model, hve parameters determine 
the hopping probabilities for the particles. In the bulk, the probability of hopping left is q times the 
probability of hopping right. The four other parameters refer to the probabilities for a particle to 
enter and exit at the left and right boundaries of the lattice. This model has attracted considerable 
attention in physics, and in particular is treated as a toy model in the study of dynamical systems 
far from thermal equilibrium. The ASEP can be solved completely in the general case (i.e. there are 
explicit formulas for the stationary probabilities) and it exhibits boundary-induced phase transitions. 
A seminal tool is the so-called Matrix Ansatz of Derrida, Evans, Hakim and Pasquier |3]. 

In the past ten years, much work has been done in combinatorics in relation to the ASEP. Most 
notably, the stationary probabilities can be interpreted as a weighted sum over certain combina¬ 
torial objects. The main starting point for the combinatorial interpretations was the ASEP with 
three parameters a, f3 and ( 7 , i.e. the partially asymmetric exclusion process (PASEP) (this is the 
same model as the ASEP, except that some parameters are set to 0). In this process, particles can 
enter the lattice on the left with probability proportional to a and exit at the right with probability 
proportional to /?, and hop right or left inside the strip with probabilities proportional to 1 and q 
respectively. In relation with the Matrix Ansatz various authors have introduced several equiva¬ 
lent combinatorial objects such as permutation tableaux, alternative tableaux, staircase tableaux, and 
tree-like tableaux, all enumerated by n! and thus in bijection with permutations. In the case where 
q = 0 (TASEP), the corresponding tableaux are enumerated by Catalan numbers. A very rich area 
of combinatorics has emerged from this topic. 

The PASEP has also been generalized to allow different species of particles, with some priority rules 
where particles can swap with particles of different species. Again, some combinatorial interpreta¬ 
tions are starting to be discovered, using an extended Matrix Ansatz. In this paper we consider 
a simple two-species PASEP with three parameters a, and q, analogous to the ordinary PASEP 
(see El HE]). In j^, the hrst author studied the case of the two-species PASEP for q = 0, and 
introduced an object called the “multi-Catalan tableau” that gives an interpretation for the steady 
state probabilities of the process for that special case. In this work, we generalize the result for all 
q with a new object called the rhombic alternative tableaux. 
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More specifically, the two-species PASEP has two species of particles, one heavy and one light. The 
“heavy” particle can enter the lattice on the left with rate a, and exit the lattice on the right with 
rate (3. Moreover, the “heavy” particle can swap places with both the “hole” and the “light” particle 
when they are adjacent, and similarly the “light” particle can swap places with the “hole”. Each of 
these possible swaps occur at rate 1 when the heavier particle is to the left of the lighter one, and at 
rate q when the heavier particle is to the right (we simplify our notation by treating the “hole” as a 
third type of “particle”). More precisely, if we denote the “heavy” particle by D, the “light” particle 
by A, and the “hole” by E, and let X and Y be any words in {D, A, E} then the transitions of this 
process are: 

AAET^iAEAT ADET^AEDT ADAy^iAADT 

<? 9 <? 


EA-^DA AD^AE 

where by A A we mean that the transition from A to A has probability n being the length 
of A (and also A). 

Notice that since only the “heavy” particle can enter or exit the lattice, the number of “light” particles 
must stay fixed. In particular, if we fix the number of “light” particles to be 0, we recover the original 
PASEP. 

Uchiyama provides a Matrix Ansatz along with matrices that satisfy the conditions, to express the 
stationary probabilities of the two-species PASEP as certain matrix products. 

Theorem 1.1 (Uchiyama [9]). Let W = {VPi,..., Wn} with Wj G {D, A, E} for 1 < i < n represent 
a state of the two-species PASEP of length n with r “light” particles. Suppose there are matrices D, 
E, and A and vectors (tc| and |r;) which satisfy the following conditions 


DE = D + E + qED DA = A-\- qAD AE = A-{- qEA 


(wlE = —{w 
a 


D\v) = -\v) 


then 


Prob(r) — ——(-wl D l(rv.=D) l(Wi=A) +E l(Wi=E) 1^^) 

^n.r 1 

’ t=\ 

where Zn,r is the coefficient of y” in • 


This result generalizes a previous Matrix Ansatz solution for the regular PASEP of Derrida et. al. in 
|1]. Inspired by Uchimaya’s Matrix Ansatz, we introduce in this paper rhombic alternative tableaux 
(RAT). They are defined in Sectionbut we state our main theorem below. 


Theorem 1.2. Let W he a state of the two-species PASEP of size n with exactly r “light” particles. 
Then the stationary probability of state W is 


Pr(lA) = ^^wt(r) 

^n,r rr 


where T ranges over the rhombic alternative tableaux corresponding to W, wi{T) is the weight of 
such a tableau, and Zn,r is the weight generating function for the set of rhombic alternative tableaux 
corresponding to the state space ofW. 
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For the rest of the paper, in Section]^ we introduce the rhombic alternative tableaux, and in Section 
[^we prove Theorem 1.2 In Section]^ we conclude with a discussion of a special case for the rhombic 
alternative tableaux when q = 0. 


2 Rhombic alternative tableaux 


The rhombic alternative tableaux (RAT) are an analog on a “triangular lattice” of the 
alternative tableaux |in) that correspond to the ordinary PASEP. By triangular lattice, 
we mean one which has as its vertices the integer points and the possible edges are 

the south edges with vertices {i,j — 1)}, west edges with vertices {(i, j), (i — 

and southwest edges with vertices {(i, j), (i — l,j — 1)} for integers i,j, as in the figure on 



the right. 


2.1 Definition of RAT 

Definition 2.1. Let IF be a word in the letters {D, A, E} with k D’s, £ E’s, and r A’s of total length 
n := k + £ + r. Define Pi to be the path obtained by reading W from left to right and drawing a 
south edge for a D, a west edge for an E, and a southwest edge for an A. (From here on, we call 
any south edge a D-edge, any west edge an E-edge, and any southwest edge an A-edge.) Define P 2 
to be the path obtained by drawing £ west edges followed by r southwest edges, followed by k south 
edges. A rhombic diagram r(lF) of type IF is a closed shape on the triangular lattice that is 
identified with the region obtained by joining the northeast and southwest endpoints of the paths 
Pi and P 2 (see Figure]^. 



Figure 1: The tiles DE, DA, and AE. Figure 2: (a) maximal and (b) minimal hexagons. 

Definition 2.2. A tiling 7~ of a rhombic diagram is a collection of open regions of the following 
three parallelogram shapes as seen in Figure the closure of which covers the diagram: 

• A parallelogram with south and west edges which we call a DE tile. 

• A parallelogram with southwest and west edges which we call an AE tile. 

• A parallelogram with south and southwest edges which we call a DA tile. 

We define the area of a tiling to be the total number of tiles it contains. 

By convention, we label the E-edges of the southeast boundary of the rhombic diagram with 1 
through £ from right to left, and the D-edges with 1 through k from top to bottom. 

Definition 2.3. A north-strip on a rhombic diagram with a tiling is a maximal strip of adjacent 
tiles of types DE or AE, where the edge of adjacency is always an E-edge. A west-strip is a maximal 
strip of adjacent tiles of types DE or DA, where the edge of adjacency is always a D-edge. The i’th 
north-strip is the north-strip whose bottom-most edge is the i’th (from right to left) E-edge on the 
boundary of the rhombic diagram. The j’th west-strip is the west-strip whose right-most edge is 
the j’th (from top to bottom) D-edge on the boundary of the rhombic diagram. Figure]^ shows an 
example of the west- and north-strips. 
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Note that the number of tiles in the i’th north-strip is the total number of D’s and A’s in the word 
W preceding the i’th E. Similarly, the number of tiles in the j’th west-strip is the total number of 
E’s and A’s in the word W following the j’th D. 


2 1 



Figure 3: Left) west-strips and (right) north-strips. 


Finally we define the main object we are working with. 

Definition 2.4. A RAT of type VE is a rhombic diagram r(VE) and an arbitrary tiling T with DE, 
DA, and AE tiles, and a filling E of T with a’s and /3’s with the following conditions: 

i. A DE tile is empty or contains an a or a /3. 

ii. A DA tile is empty or contains a /3. 

iii. An AE tile is empty or contains an a. 

iv. Any tile above and in the same north-strip as an a must be empty. 

V. Any tile to the left and in the same west-strip as a /3 must be empty. 

We define fi(fT, T) to be the set of fillings of tiling T of the rhombic diagram r(VE). In other words, 
F G fi(VE, T) means E is a filling of type W of the tiling E. 

Definition 2.5. A north line is a line drawn through each north-strip containing an a, starting 
at the tile directly above that a. A west line is a line drawn through each west-strip containing a 



Figure 4: An example of a RAT of size (9,3,4) 
with type DAADDEDAE and weight 


Figure 5: A complete representation of a RAT 
that is equivalent to the example on the left. 
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13, starting at the tile directly left of that /3. An example of the RAT with the north- and west lines 
is shown in Figure 


In terms of the north- and west lines, we rewrite the 
conditions (iv) and (v) of Definition 2.4 by (equivalently) 
requiring that any tile that contains a north line or a 
west line must be empty. 


Definition 2.6. The size of a RAT of type W is 
{n,r,k), where k is the number of D’s in W, r is the 
number of A’s in W, and n is the total number of let¬ 
ters in W. We can also call this the size of a filling F 
of type W. We can also refer to the size of a tableau as 
simply (n, r), where we do not keep track of the number 
of D’s. 



Definition 2.7. To compute the weight wt(F) of a Figure 6: The minimal tiling of a rhombic 
filling F, first a g is placed in every empty tile that does diagram r(D AAEEDEAE) defined by 
not contain a north line or a west line. Next, wt{F) is the paths Pi and P 2 . 
the product of all the symbols inside F times a^f3^, for 
F a filling of size (k + i + r,r,k). 

We will prove in Proposition |2.8| that the sum of the weights of all fillings of r(VF) does not depend 
on the tiling F. 


2.2 Independence of tilings and definition of weight(kF) 


Proposition 2.8. Let W be a word in {D, A, E}. Let Ti and T 2 represent two different tilings of a 
rhombic diagram r(VF) with DE, DA, and AE tiles. Then 

Feaiw,Ti) f’&&{w,T2) 


Definition 2.9. Consider a hexagon with vertices {(i,j), {i,j — 1), (i — l,j — 2), (z — 2,j — 2), (i — 
2, j — 1), (i — 1, j)} for some integers i, j that is tiled with a DE-, a DA-, and an AE tile. A minimal 
hexagon is when the tiles have the configuration of Figure]^ (b), and a maximal hexagon is when 
the tiles have the configuration of Figure]^ (a). 

Definition 2.10. Let IF be a word in {D,A, E}. We define the minimal tiling of r(lF) to be 
the tiling that that does not contain an instance of a maximal hexagon, such as the example in 
Figure We refer to such a tiling by Fnin can be constructed by placing tiles from Pi 

inwards, and placing an AE tile with first priority whenever possible. (Similarly, a maximal tiling 
is one that that does not contain an instance of a minimal hexagon, and is referred to by Fnax- The 
maximal tiling can be constructed by placing tiles from Pi inwards with priority given to the DA 
tiles whenever possible.) 


Definition 2.11. A flip is an involution that switches between a maximal hexagon and a minimal 
hexagon, and is the particular rotation of tiles that is shown in Figure 


^In the remark following the proof of Lemma 


2.12 


we show that Tmin is the unique minimal tiling. 
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The lemma below contains a generally known result, no¬ 
tably in the case of a plane partition. 

Lemma 2.12. Let r(lT) he a rhombic diagram of type W. 

For any two tilings T and S ofT{W), T can he obtained 
from S by some series of flips. 

Remark. To make the paper self-contained, we will give here a proof that tilings are in bijection with 
configurations of non-crossing paths, a well-known fact. Our proof dehnes a classical construction, 
in particular in the case of a plane partition where the bijection will give a configuration of paths 
related to a binomial determinant (by the LGV lemma), which can expressed by a simple formula 
giving the well-known MacMahon formula for plane partitions (or 3D Ferrers diagrams) within a box 
of size (a, b, c). We note that the case of a plane partition within a box of size {k, r, i) is equivalent 
to the tilings of r(iy) where W = D^A^E^. 

Definition 2.13. A non-crossing configuration {r} is a collection of lattice paths {ri,r 2 ,...} 
where no two paths share a vertex (and hence are non-crossing). 



E E 


Figure 7: A flip from a maximal (left) 
to a minimal hexagon (right). 


Proof of 2.1S\ We prove the lemma by constructing a bijection from certain non-crossing conhg- 
urations on r(VF) to tilings of r(lF). We define a height function ht({T}) on the non-crossing 
configurations. We dehne the minimal non-crossing conhguration {T^rain to be the one with height 
0. We also dehne a local move called a slide on {r} which we show corresponds to a hip on T (sim¬ 
ilarly, a reverse slide corresponds to a hip in the other direction). We show that a slide diminishes 
the height of a non-crossing conhguration by 1, and that any non-crossing conhguration with height 
greater than 0 admits a slide. Thus it is possible to obtain {rjmm from any {r} by some set of 
slides, and so it is possible to obtain any {cr} from any {r} by some set of slides and reverse slides. 
This translates to the desired result due to the bijection. 



Figure 8: (Left) a tiling T with west-strips indicated, bijection to the non-crossing conhguration 
{r} with ht({r}) = 2. (Right) a slide performed at the indicated free lattice point to obtain the 
non-crossing conhguration {r'} with ht({r'}) = 1, which is in bijection with a different tiling E' 
(also obtained from T by performing a hip at the northeast-most hexagon). 


Bijection from tilings on r(VF) to non-crossing configurations on r(VF). Recall that r(kF) is on a 
lattice with integer points {(^,j)}, and a tiling on T{W) consists of tiles DE, DA, and AE whose 
corners are on those lattice points. We dehne a bijection from a tiling 7~ to a non-crossing conhgu¬ 
ration {r} on the lattice contained within r(kF). Let {r} = {rj,... ,rfc}, where each r* is a lattice 
path consisting of west and southwest steps {j, t) — )■ (j — 1, i) and (j, i) — )■ {j — 1,£— 1 ) respectively. 
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that starts at e* and ends at vt for some sets {cj} and {vi}. We let be the north endpoint of 
the i’th D-edge on the southeast boundary of r(W) (from top to bottom), and let Vi be the north 
endpoints of the z’th D-edge on the northwest boundary of r(W) (from top to bottom). 

To obtain {r} from T, let r* be the path that coincides with the northwest boundary of the i’th 
west-strip (from top to bottom). Since the west-strips do not cross each other, it is clear that {r} 
is well defined in this way (see Figure]^. 

Define a free lattice point to be a lattice point that is not part of any path in {r}, and does not lie 
on the southeast boundary of r(VF). To obtain T from {r}, we do the following: for each r^, place 
a DE tile directly below and adjacent to each west step of r,, and place a DA tile directly below and 
adjacent to each southwest step, so that the northwest boundaries of the DE and DA tiles coincide 
with the steps of r*. Thus Tj corresponds to a west-strip. Now, for every free lattice point that is 
not on the southeast boundary of r(bF), place an AE tile so that its northwest corner coincides 
with that lattice point. We claim that we obtain in this way a valid tiling of r(VF). To check this, 
we must simply verify that the construction above results in no tiles overlapping, and a complete 
covering of the shape. 

Height of {r}. For each i, define the minimal path mi to be the one starting at and taking a 
maximal possible number of west steps followed by a maximal possible number of southwest steps 
to Vi- Define the height of Ti G {r} (i.e. ht(Tj)) to be the area between r* and m^ (i.e. the number of 
free lattice points strictly northwest of r* and weakly southeast of mi). We define the height of {r} 
to be 

ht({r}) = ^ ht(ri). 

It is clear that there is a unique {r} of height 0, by letting each r* be mi. We call this {rmm}- 

A slide on {r}. Let p be a free lattice point {i,j) such that (i -|- l,j), {i,j — 1), and (i — l,j — 1) 
are all not free lattice points. Then, those lattice points must necessarily belong to the same 
path in the non-crossing configuration {r}, say r*. A slide on {r} at the location of p means 
exchanging the steps southwest and west for the steps west and southwest in Ti, and thereby passing 
through p and creating a new free lattice point below Ti, as in Figure More precisely, the steps 
{i + l,j) —t {i,j — 1) —)> {i—l,j — 1) in Ti are exchanged for {i + l,j) —>■ {i,j) —)• — — 1) to make 

Tjb Clearly does not cross {ri,..., rj_i, Tj+i,...}, and so the new collection of paths {r'} formed 
by replacing r* with is also a non-crossing configuration. Furthermore, since ht(TjQ = ht(rj) — 1, 
we have ht({r'}) = ht({r}) — 1 . 

We have already established that a southwest step of r* in {r} corresponds to a DA box in the west- 
strip i, and a west step of Ti corresponds to a DE box. A free lattice point necessarily corresponds to 
an AE tile, since all DA and DE tiles must be part of some west-strips. From Figure]^ it is easy to 
see how a slide corresponds to a flip from a maximal hexagon to a minimal hexagon, and a reverse 
slide (the inverse operation) corresponds to the reverse flip. 

Notice that if no such free lattice point p exists such that its three neighbors east, south, and 
southwest are all not free lattice points, this implies r* = m, for each i. Then {r} is the minimal 
non-crossing configuration, as in Figure]^ Thus {r} admits no slides if and only if it equals {Tmin}- 

Now we complete the proof. Let {r} be some non-crossing configuration with ht({r}) = A; > 0. 
Then by the above there is at least one free lattice point that admits a slide. After performing a 
slide at that location, we obtain a new non-crossing configuration {r^} with ht({r'}) = ht({r}) — 1 . 
Recursively, this implies that by applying some series of slides, we can get from {r} to a non-crossing 
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Figure 10: A more explicit set of tiles that is in simple correspondence with the RAT fillings. 


configuration with height 0. However, {rmin} is the unique such non-crossing configuration, so we 
have shown that we can get from {r} to {Tmin} with some set of slides. We can now equivalently 
define ht({r}) as the minimal number of slides required to get from {r} to {Tmin}- 

Let {(t} be a different set of non-crossing configurations. There is similarly a set of slides to get from 
{cj} to {Tmin}- Thus we Can get from {r} to {fi} by a series of slides, by first applying the slides 
to get from {r} to {Tmin}, and then by applying slides in reverse to get from {Tmin} to {cr}. Let 
the tiling T correspond to the non-crossing configuration {r}, and the tiling S to the non-crossing 
configuration {u}. Since the slides on the paths correspond to flips on the tilings, we obtain that 
one can get from T to 5 with a series of flips, as desired. □ 

Remark. It is easy to check that {Tmin} corresponds to the minimal tiling Tmin of r(lT) as in Figure 
since Tmin is a tiling that admits zero flips from a maximal hexagon to a minimal hexagon, 
which corresponds to a non-crossing configuration admitting zero slides and having height 0. Since 
{Tmin} IS the unique non-crossing configuration of height 0, Tmin must be the unique minimal tiling 
according to our definition of “minimal”. (The maximal tiling Tmax is also unique by a similar 
argument.) Furthermore, the minimal number of flips required to get from T to Tmin is commonly 
referred to as the height of a tiling T. 

For the proof of Lemma [2.8[ we introduce a more 
explicit set of tiles, where tiles can now contain 
a, /?, q, a north line, a west line, or both a north 
line and a west line, as in Figure We can 
now describe a complete covering of S with the 
following. 

Definition 2.14. Two adjacent tiles are com¬ 
patible if: 

1 . a tile has a north line through it if and 
only if its south E-edge is adjacent to a 
tile containing an a or a north line, and 

2 . a tile has a west line through it if and only 
if its east D-edge is adjacent to a tile con¬ 
taining a /3 or a west line. 



Figure 9: Tmin and 


Proof of 2.8 For any rhombic diagram r(VF), any two tilings 71 and T can be obtained from each 
other by some series of flips by Lemma [2.12[ Thus it is sufficient to show that if 71 and T differ by 
a single flip, then there is a weight-preserving bijection between fi(LF, 71) and fi(fT, Ti). Let this flip 
occur at a certain “special hexagon” (f)i in 71 and 1)2 m 71). Without loss of generality, let l)i be of 
minimal type as on the left of Figure and let 1)2 be of maximal type as on the right of Figure 
The rest of the tiles are identical in 71 and 71. 


We define the bijection from F G fi(lT, 71) to some F' G fi(fT, 71) with an involution cf). To begin. 
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(j) sends every tile including its contents in 71\fii to its identical copy in 72\1)2- Then, (p sends the 
tiles and contents of f)i to a rearrangement of those tiles according to the cases shown in Figure pT| 
It is easy to see that this map preserves the weights of the hllings, since the quantities of /3s, a’s, 
and q's are preserved for each case. 


We claim that the map p also preserves the compatibility of the tiles, as dehned in Dehnition 2.14 
It is easy to check that in each possible case of f}i, the tile adjacent to the south E-edge contains 
an a or a north line if and only if the tile adjacent to the south E-edge of contains an a or 

a north line. Similarly, the tile adjacent to the east D-edge of f)i contains a /3 or a west line if and 
only if the tile adjacent to the east D-edge of p{\]i) contains a ,0 or a west line. 


Thus p indeed gives a weight-preserving bijection from fi(iy, 71) to fi(IT, Ti), and so the lemma 
follows. □ 


Thus we are able to make the following dehnition. 

Definition 2.15. Let LF be a word in {D, A,E}, and let T be an arbitrary tiling of r(LF). Then 
the weight of a word W is 

weight(IF) = ^ wt(F). 

Fefi(iv,r) 



Figure 11: The involution p from each possible hlling of a minimal hexagon (left) to a maximal 
hexagon (right). The arrows imply compatibility requirements. 


In fact, we can dehne equivalence classes of rhombic alternative tableaux with the following dehni- 
tions. 


Definition 2.16. A weight-preserving flip on a RAT with tiling T is the transformation p (or 
the inverse of the p) given by Figure 11 on some hexagon f) of T and the symbols contained in it, 
while preserving the hlling of T\f). 


Definition 2.17. Let LF be a word in {D,A, E}, and let Ti G h(LF, 71) and T 2 G h(LF, 72 ) be 
rhombic alternative tableaux of type W for arbitrary tilings 71 and Ti of r(VF). Then Ti and T 2 
are equivalent if and only if T 2 can be obtained from Ti by a series of weight-preserving hips. 
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2.3 Enumerative considerations 


In this subsection, we show that the number of RAT equivalence classes of size (n, r) is (”) 

Definition 2.18. Define 

Zn,r(a,l3,q) = ^weight(iy) 

w 

for W ranging over all words in {D, A, E}” with r A’s. 

By convention, let Zn,n{c(, P, q) = 1- 

Theorem 2.19. 

n—1 

Y\{o( + (3 + iaP). ( 2 . 1 ) 

i=r 

Proof. Let Zn^r,k{ci, f3,q) be the weight generating function for the RAT (with the maximal tiling) 
with exactly k west-strips that do not contain a (3. 

We also define 

Zn,r{x) = ^ Zn,r,k{ot, / 3 , l)x^ ( 2 . 2 ) 

k>0 

with Zn,r{o:,f3, 1) = Zn,r{^)- We claim that 

n— 1 

+ (3 + iaj3). (2-3) 

i=r 



Zn,ria,l3, 1 ) = 


n 


We will prove Equation (2.3) and hence Equation (2.1) by induction on n in terms of the more 
refined Zn^r,k^- 

First, when n = 1, Zi^o{x) = xa + (3 and Zip = 1 by convention. 


Now, we suppose that Equation (2.3) holds for any N < n and any r < n — 1. (Again, by convention 
Zn,n = 1 for all n.) We will show that the formula holds as well for N = n -|- 1, for all r < n. 


We begin by observing that Equation (2.3) satisfies 


Zn+l,r{x) = (xa + P + ra/3)Zn,r{x + /3) + Zn,r-l{x + f3). 


(2.4) 


We now construct a recursion for Zn+i^r,k in terms of the functions {Z^y^k'} by keeping track of the 
terms after the addition of a new letter D, A, or E to the end of a word W G {D, A,E}"'. For the 
following, we denote by T' a tableau in ^^iW,Tmax{W^)) (for TmaxiW) the maximal tiling of r(lT)), 
and by r a tableau in fi{Wx^Trnax{Wx)) for x G {D,A, E}. We consider all possible cases for W 
and corresponding T' such that the resulting T G ^{Wx,Tmax{Wx)) has size (n -|- l,r) and exactly 
k west-strips that do not contain a (3. 


1. We add a D to the end of a PASEP word W of length n with r A’s. On the tableau 
level, this corresponds simply to the addition of a D-edge to the southwest end of each 
T' G ^{W,%nax{W)) as in Figure 12 (a). Since the hlling of T is the same as that of T', 


Z, 


n+l,r,A: 


= aZ, 


n,r,/c—1 • 


(2.5) 


This contributes to wt(T)|(g=i). 
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Figure 12: Adding a (a) D, (b) A, or (c) E to the end of W. 


2. We add an A to the end of a PASEP word W of length n with r — 1 A’s. On the tableau level, 
this corresponds to the addition of a vertical column of some DA tiles to the left boundary of 


T’ to form a tableau with the maximal tiling of T{W A) as in Figure 12 (b). Suppose T' has 


i > k west-strips that do not contain a /?. Then to obtain T with exactly k west-strips that 
do not contain /3, the ^ — k DA tiles that do contain a (3 can be chosen in ways, with the 
other k DA tiles containing q. Therefore, we obtain 


Pi f ) 


E 

t>k 


P 


e-k 


y 

^n,r— 


!,£■ 


( 2 . 6 ) 


This contributes to A))^^{'^)\{q=iy 

We add an E to the end of a PASEP word W of length n with r A’s. On the tableau level, 
this corresponds to the addition of a vertical column of some DE tiles followed by a strip of r 
AE tiles to the left boundary of T' to form a tableau with the maximal tiling of r(VFE) as in 
Figure 12 (c). Suppose T' has i > k west-strips that do not contain a p. We have two cases. 


(1) For the first case, there is no o. in the newly added DE tiles. Then to obtain T with 
exactly k west-strips that do not contain ,0, the i — k DE tiles that do contain a P can be 
chosen in (^) ways, with the other k DE tiles containing q. Following this, the AE tiles 
can either contain all q’s, or some consecutive string of q’s followed by an a. 

(2) For the second case, there is an a in the newly added DE tiles, with some l!. — k<u<i — l 

“free” DE tiles below it. (A “free” tile means the tile lies in a west-strip that does not 
contain a p.) Then, to obtain T with exactly k west-strips that do not contain P, the 
i — k DE tiles that do contain a /3 can be chosen in ways, with the other u — {£ — k) 

DE tiles that lie below the a containing q, and the rest of the tiles empty. 

Combining the above two cases, we obtain 


^n+l,r,k{.^i Pi^} P 


E 

e>k 


P^ ^{ra + l)Zj^r,i + 


i-i 

E 

u=£—k 


This 


contributes to ^T£&{WE,Tmax(WE)) wt(r)|(,=i). 


u 

£-k 


aP' 


i—k y 
Zj. 




(2.7) 
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Combining Equations (2.5), (2.6), and (2.7) and summing over k, we obtain 


Zn+l,ri^) = ^ ( XaZ^,r,k-l + ^ (^^P^-’^x'^Zn,r-l,i + /? (^^^^-’^X^Zn,rdra + 1 ) 

fc>0 \ £>k ^ ^ V / 

£-1 

+/^E E 


£>k 


i>k u=£—k 




xaZn,r{x) + Zn,r-i{x + /?) + l3{ra + l)Zn,r{x + /3) + a/3 eeu ! i )< 3 '-‘-‘ 2 „.,, 

k>0£>k ^ ^ 

XaZn^rix) + Zn^r-\{x + (3) + {raj3 + (3 + Xa)Zn^r{x + (3) — xa E 


A:>0 


which simplihes to Equation (2.4). Since satishes the desired recursion, we thus obtain that 

Equation (2.3) indeed holds for N = n + 1, and so our proof is complete. □ 


Remark. It would be interesting to obtain a bijective proof for Equation (2.1). 


3 Steady state probabilities of the two-species PASEP 


Our main result is the following. 


Theorem 3.1. Let W be a word in {D, A, E}” that represents a state of the two-species PASEP 
with exactly r “light" particles. The stationary probability of state W is 


Pr(IT) = —^—weight(IT), 

^n..r 


where weight(VE) is defined in Definition 2.15 




(3.1) 



D 


Figure 13: 


All hllings for the minimal tiling of a rhombic diagram of type DAE. 


Example. All seven hllings of the minimal tiling of a rhombic diagram corresponding to the word 
DAE are shown in Figure From the sum of the weights of these hllings, we obtain 

Pr(DAE) = —— + aq^ + aq + (3q‘^ + Pq + Oi/3 + a/dq) . 

Z3,l 

To facilitate our proof, we provide a more hexible Matrix Ansatz that generalizes Theorem o with 
the same argument as in an analogous proof for the ordinary PASEP of Corteel and Williams [3l 
Theorem 5.2]. For a word W G {D, A,E}’^ with r A’s, we dehne unnormalized weights f{W) which 
satisfy 

Pr(iy) = f{W)IZn,r 

where Zn,r = 'T^w' /(^O where the sum is over all words W' of length n and with r A’s. 
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Theorem 3.2. Let X be a constant. Let (t(;| and |i)) be row and column vectors with (u)||u) = 1. Let 
D, E, and A be matrices such that for any words X and Y in {D, A, E}, the following conditions 
are satisfied: 

(I) {w\X{DE - qED)Y\v) = X{w\X{D + E)Y\v), 

(II) {w\X{DA - qAD)Y\v) = X{w\XAY\v), 

(III) {w\X{AE - qEA)Y\v) = X{w\XAY\v), 

(IV) /3{w\XD\v) = A(u;|A|?;), 

(V) a{w\EY\v) = A(u;|y|z;). 

Let W = {Wi,..., Wn} for Wi G {D, A, E} for 1 < i < n. Then for any state W of the two-species 
PASEP of length n with r “light" particles, 

1 

^ {w\A^\v) +Al^Wi=A) +E1(w,=e)\v). 


Proof. The proof of Theorem 3.2 follows exactly that of |3l Theorem 5.2], Note that the above 
implies that 

{w\{D + yA-\- E)'^\v) 


Zn,r = [yl- 




□ 


3.1 Proof of the main theorem 


Proof of 3.1 The Matrix Ansatz of Theorem |3.2| implies that the steady state probabilities for 
the two-species PASEP satisfy certain recurrences (that in turn determine all probabilities). The 
strategy of our proof is to show that the weight generating function for RAT of fixed type satisfies 
the same recurrences. Specifically, we use these recurrences with the constant A = ajS, to show by 
induction that for W a word in {D, A, E}” with r A’s, 


f{W) = weight(lE). 


(3.2) 


We do the induction on the area M of our tableaux. 


By Theorem 3.2, f{W) = 


VE as a product 


1 


{w\A^\v) 
of the matrices D 


(r(;|lE|u) for W G {D, A, E}” a word with r A’s, where we interpret 


E, 


and A in the order those letters occur in IE. To start, if the 
area of r(lE) is 0, then necessarily W = E^A'’D^ for some I and k. Eor every such case, there is 
a single tiling of r(lE) and a single RAT on that tiling, and both of these a re tr ivial (i.e. there are 
zero tiles to be filled.) We obtain weight(E^ A'" D^) = 


From Theorem 


base case indeed satisfies Equation (3.1) 


3.2 


it is clear that the 


Now suppose that any word of area M < m satisfies Equation (3.1). Let VE be a word of length n 
with r A’s, such that r(lE) has area m. Outside of the base case, we assume that at least one of 
the following must occur: 

i. VE contains an instance of “DE”. 


ii. VE contains an instance of “DA”. 
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iii. W contains an instance of “AE”. 

Based on the occurrence of one of the above, we will express weight(Ty) in terms of the weight of 
some other words whose rhombic diagrams have areas smaller than m. Throughout the following, 
we let X and Y represent some arbitrary words in {D, A, E}, and we let T represent a RAT of type 
W. 


(i.) W contains an instance of “DE”. We write W = ADET. We can choose an arbitrary tiling 7~ 
of r(VE), since any such tiling will contain a DE tile adjacent to the chosen DE edges. We call this 
DE corner tile the chosen corner. Let T G fi(W, T). The chosen corner of T must contain either an 
a, a /3, or a q, so we can decompose the possible fillings of T into three cases. 


If the chosen corner contains an a, then all the tiles above it in the same north-strip are empty, and 
so its entire north-strip has no effect on the rest of the tableau. Thus such T can be mapped to a 
filling of a smaller RAT on tiling E' with that north-strip removed, which would have type XDY 


(similar to the example in Figure 14 (a)). It is easy to check that this operation results in a valid 


tableau, since any two symbols in the same west-strip of T remain in the same relative position in a 
west-strip of T'. (And similarly for the north-strips, save for the one that was removed). This map 
gives a bijection between tableaux of type A D E T on tiling T with an a in the chosen DE corner 
and tableaux of type ADA on tiling T' ■ The removed column with the a in its bottom-most box 
has total weight a/ 


Similarly, if the chosen corner contains a /3, then the tiles to its left in the same west-strip must be 
empty, and so its entire west-strip has no effect on the rest of the tableau. Hence such T can be 
mapped to a smaller RAT on tiling T" with that west-strip removed, which would have type A E A 


as in Figure 14 (b). This map gives a bijection between tableaux of type ADEA on tiling T with 
a /3 in the chosen DE corner and tableaux of type A E A on tiling E"■ The removed west-strip with 
the (5 in its right-most tile also has total weight af3. 


Finally, if the chosen corner contains a g, then this tile has no effect on the rest of the tableau. Hence 
such E can be mapped to a RAT of area m — 1 on tiling E'" with that DE corner tile removed, 
which would have type A ED A (similar to the example in Figure [^(c)). This map gives a bijection 
between tableaux of type ADEA on tiling E with a g in the chosen DE corner and tableaux of 
type A ED A on tiling E'"■ The removed tile with the g has total weight g. 


Consequently, we have the sum of the weights of the fillings: 


weight(ADEA) = weight(AD A) • a/3 -|- weight(AEA) • a/3 -|- gweight(AED A). 

By the induction hypothesis, since the areas of r(ADA), r(AEA), and r(AEDA) are all 
strictly smaller than m, we have weight(ADA) = /(ADA), weight(AEA) = /(AEA), and 
weight(AED A) = F(AEDA). Thus we obtain 


(a/3)(rc|A(D + E)Y\v) + q{w\XEDY\v) = {w\XDEY\v) = {w\W\v). 


Hence by Theorem 3.2 with A = af5, it follows that W satisfies Equation (3.2). 


(a.) W contains an instance of “AE”. We write W = AAEA. We choose tiling E of r(IA) such 
that there is an AE tile adjacent to the chosen AE edges (and we allow the rest of the tiling to be 


^In the total weight of a column, we include the weight of the bottom-most edge, which is a component of the 
southeast boundary of T. When the column removed is an E-column, the weight of the boundary component is /3, so 
the total weight of the column with an a at the bottom is otfi. Similar reasoning is used in the other cases. 
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Figure 14: (a) X AEY ^ XAY, (b) XBEY ^ XDY, and (c) XDAF qXADY. 


arbitrary). We call this AE corner tile the chosen corner. Let T G fi(VF, 7~). The chosen corner of 
T must contain either an a or a 

If the chosen corner contains an a, then the tiles above it in the same north-strip must be empty, 
and so its entire north-strip has no effect on the rest of the tableau. Hence such T can be mapped 
to a smaller RAT on tiling T’ with that north-strip removed, which would have type X AY, as in 


Figure 14 (a). This map gives a bijection between tableaux of type X AEY on tiling 7~ with an a 
in the chosen AE corner and tableaux of type A AT on tiling Y'. The removed north-strip with 
the a in its bottom-most tile has total weight a 15. 

On the other hand, if the chosen corner contains a q, then this tile has no effect on the rest of the 
tableau. Hence such T can be mapped to a RAT of area m — 1 on tiling Y" with that AE corner 


tile removed, which would have type XEAY (similar to the example in Figure 14 (c)). This map 


gives a bijection between tableaux of type X AEY on tiling Y with a g in the chosen AE corner 
and tableaux of type XEAY on tiling Y". The removed tile with the q has total weight q. Thus 
we obtain the sum of the weights of the fillings: 

weight(A AET) = weight(A AT) • a/5 -|- gweight(AEAy). 

Similar reasoning to the DE case completes the argument. 

(Hi.) W contains an instance of “DA”. We write W = XE> AY. We choose tiling Y of r(VF) such 
that there is a DA tile adjacent to the chosen DA edges (and we allow the rest of the tiling to be 
arbitrary). We call this DA corner tile the chosen corner. Let T G fi(W,T). The chosen corner of 
T must contain either a /? or a g. 

If the chosen corner contains a /3, then the tiles to its left in the same west-strip must be empty, 
and so its entire west-strip has no effect on the rest of the tableau. Hence such T can be mapped 
to a smaller RAT on tiling Y' with that west-strip removed, which would have type X AY (similar 


to the example in Figure 14 (b)). This map gives a bijection between tableaux of type AD AT on 
tiling Y with a (5 in the chosen DA corner and tableaux of type A AT on tiling Y'. The removed 
west-strip with the (5 in its right-most tile has total weight al5. 

On the other hand, if the chosen corner contains a q, then this tile has no effect on the rest of the 
tableau. Hence such T can be mapped to a RAT of area m — 1 on tiling Y" with that DA corner tile 


removed, which would have type A ADA, as in Figure 14 (c). This map gives a bijection between 


tableaux of type A D A A on tiling Y with a q in the chosen DA corner and tableaux of type A A D A 
on tiling Y". The removed tile with the q has total weight q. Thus we obtain the sum of the weights 
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of the fillings: 


weight(XDAy) = weight(XAY) • aj3 + gweight(AADI"). 

Similar reasoning to the DE case completes the argument. 

From the above cases, we obtain that for any M, any word W with r(bE) of area M satisfies 
Equation ( |3.2[ ), which is the desired result. □ 

An independent proof of Theorem |3.1| can be obtained by constructing a Markov chain on the rhombic 
alternative tableaux that projects to the two-species PASEP, and is featured in a forthcoming paper 
by the first author. 


3.2 Complements: the two-species PASEP algebra and the cellular ansatz. 


It would be possible to give another proof of Theorem 3.1 in the spirit of the general theory called 
cellular ansatz, introduced and developed by the second author in m- 


In the simple case r = 0 of the PASEP, the Matrix Ansatz defines an algebra with generators D 
and E, with relation DE = g'ED-|-E-|-D. In this algebra, any word W with letters D and E can 
be written in a unique way as a sum of monomials g* D* E-^. 


The proof relies on a planarization of the rewriting rules DEi—)-( 7 ED, DEi— t-E, and DE i—?■ D 
(see an example on slides 25-50 of Chapter 3a of El). In this context, alternative tableaux appear 
naturally. We obtain an identity expressing the word W as where the sum is over 

alternative tableaux T, and wt(T) is a certain monomial of the form g*D*E-^, where i, j, and t are 
defined from the tableau T (see slide 59 of Chapter 3a of m)- By applying the Matrix Ansatz 
for the PASEP, we get immediately the interpretation of the stationary probabilities in terms of 
alternative tableaux (slide 60 of Chapter 3a of |11)1. 

The general theory of the cellular ansatz works with some family of quadratic algebras Q, having 
two families of generators, with some commutations relations. Any word W in those generators 
can be expressed as a sum of monomials over generalized tableaux called complete Q-tableaux, in 
bijection with Q-tableaux (see Chapter 6 a, slides 41-46, 54-56 of m)- 

The Matrix Ansatz for the two-species PASEP defines an algebra with three generators D, E, A and 
three commutation relations DE = gED-|-E-|-D, DA = ( 7 AD-I-A, and AE = gEA-|-A. This 
quadratic algebra does not quite fit in the general cellular ansatz theory of Chapter 6 a of HH, but 
the theory can be extended to such an algebra, by replacing the quadratic lattice by a tiling T of 
the diagram r(IT). The corresponding Q-tableaux are the RAT and, in a similar way, one can prove 
that any word W in letters {D,A, E} can be expressed in a unique way as a sum of monomials 
g* D* A™ EC 


More precisely we have the following identity 

W= ^ g‘E*A™D^ 

Fe{i(w,r) 

where i (respectively j) is the number of north-strips (respectively west-strips) of F not containing 
an a (respectively /3), and t is the number of cells weighted q as in the definition of wt(F) in Section 

is equal to the monomial 


Note that the weight wt(T) defined in Definition 


4.2 


where n is the length of W and r is the number of A’s it contains. 


Applying to the above the two-species Matrix Ansatz, we obtain immediately Theorem 3.1 
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4 Bijection from RAT with g = 0 to “multi-Catalan tableaux” 


The two-species PASEP for the case q = 0 was studied in a previous paper [7], and certain tableaux 
called “multi-Catalan tableaux” were found to give a combinatorial interpretation for the steady 
state probabilities of this process. A multi-Catalan tableau is a tableau of the same flavor as the 
alternative tableaux, except with some additional conditions on its rows and columns, which can 
be labeled as a D-row or A-row and as an E-column or A-column respectively. More precisely, the 
multi-Catalan tableaux are defined as follows. 

Definition 4.1. A multi-Catalan tableau of size {n,r,k) is a Young diagram Y = Y{T) with at 
least r inner corners (i.e. consecutive west- and south-edges on the southwest boundary), that is 
justified to the northeast and contained in a rectangle of size [k r) x {n — k). Y is identihed with 
the lattice path L = L{T) that takes the steps south and west and follows the southeast boundary 
of Y. In addition, we have the following: 

• Each boundary edge of L is labelled with a D, E, or A such that exactly r inner corners have 
both edges labeled with A’s, and the remaining west edges have the label E, and the remaining 
south edges have the label D. 

• An E-column (respectively A-column) is a column with an E (respectively A) labeling its 
bottom-most edge. 

• A D-row (respectively A-row) is a row with a D (respectively A) labeling its right-most edge. 

• A DE box is a box in a D-row and an E-column. (The DA, AE, A A boxes are defined 
correspondingly.) 

Finally, we fill T with a’s and /3’s according to the following rules: 

i. A box in the same row and left of a /3 must be empty. 

ii. A box in the same column and above of an a must be empty. 

iii. A DE box that is not forced to be empty must contain an a or a /3. 

iv. A DA box that is not forced to be empty must contain a /?. 

V. An AE box that is not forced to be empty must contain an a. 


Remark. It is easy to check that any box of the multi-Catalan tableau that has some A-row below 
it or some A-column to its left must be empty. 

In Dehnition |4.1[ when we refer to the symbol that a box “sees” to its right or below, we mean the 
first symbol encountered in the same row or column, respectively. For example, in the third tableau 
of Figure [T^ j3 is the first symbol that the (3 in the row labeled “2” sees below it. 

Definition 4.2. We call the weight wtMCTiT) of a multi-Catalan tableau T the product of all the 
a’s, /3’s it contains times where (n, r, k) is the size of the tableau. 

Definition 4.3. The type of the tableau T is the word in {D, A, E} that is obtained by reading the 
southeast boundary of L{T) from northeast to southwest, and reading a D for each D-labeled south 
edge, an E for each E-labeled west edge, and an A for each “inner corner” with a pair of A-labelled 
edges. 
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Figure 15: Bijection from a RAT of type DAADDEDAE with the minimal tiling Tmin to the “un¬ 
raveled” tableau to a multi-Catalan tableau of the same type. Here the north-strips of the RAT 
correspond to the E-columns of the multi-Catalan tableau, and the west-strips of the RAT correspond 
to the D-rows of the multi-Catalan tableau. The grey columns of the multi-Catalan tableau are the 
A-columns that were taken from the left of the “unraveled” tableau, and placed in the appropriate 
location in between the E-columns of the multi-Catalan tableau. 


Let MCT(LF) be the set of multi-Catalan tableaux of type W, and MCT(n, r) be the set of multi- 
Catalan tableaux of size {n,r). Then the following result holds for the multi-Catalan tableaux. 

Theorem 4.4 ([7]). Let W be a word in {D, A, E}"" that represents a state of the two-species PASEP 
for q = 0 with exactly r “light” particles. Let 

^n,r = ^ wt mct(T) 

TeMCT(n,r) 

Then the stationary probability of state W is 

Pr(ir) = ^ 

TeMCT(VE) 


We now provide the connection from RAT with g = 0 to the multi-Catalan tableaux. 

Proposition 4.5. Let W be a word in {D, A, E} that represents a state of the two-species PASEP. 
The multi-Catalan tableaux of type W are in bijection with the rhombic alternative tableaux with 
q = 0 onT{W) with some arbitrary fixed tiling. 


Proof. First, it suffices to prove the Proposition for any one fixed tiling, so we will do so for the 
minimal tiling Tmin of Oui' bijection follows the diagram in Figure 15 

labeling the north-strips from right to left, we perform the following: 


In particular, after 


1. “Unravel” Tmin as in Figure to obtain an “unraveled” tableau with its rows labeled D or A 
corresponding to whether the tiles in those rows came from west-strips or not (and similarly, its 
columns labeled E or A corresponding to whether the tiles in those rows came from north-strips 
or not.) Keep all a’s and /3’s in the same relative locations after unraveling. 
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2. Label the A-columns with Ai, A 2 ,... in the unraveled tableau from right to left. Label the 
E-columns in the unraveled tableau with Ei,E 2 ,. ■ ■ from right to left. Add boxes with X’s 
through them (i.e. place-holders) to correspond to squares that are both in an A-row and an 
A-column. 

3. For i = 1,2,..., take the Aj column and place it between the Ej and Ej^i column where if 
the word W is read from right to left, the Lth A is between the j’th and j + I’st E’s. For each 
row, if there is a /3 in the Aj column, move that /3 to the right-most possible empty box in the 
multi-Catalan tableau. 


We claim that the above operations result in a multi-Catalan tableau, as seen in Figure We keep 
track of the following quantities: 

• The length of column Ej (respectively Aj) is the number of D’s and A’s in W preceding the 
j’th E (respectively A). 

• The Lth D-row is composed of e* DE boxes followed by Oj DA boxes (from right to left), where 
ej is the number of E’s and Oj is the number of A’s following the f’th D in W. 


First, we show that the map we defined results in a partition with its southeast boundary edges 
labeled D, A, and E as needed according to Definition 4.1 In particular, we note that every A- 


column in the unraveled tableau has an A-row immediately below its last box since the number of 
DA boxes in consecutive D-rows is equal. This implies that whenever an A-column is inserted into 
the multi-Catalan tableau, the west A-edge on the boundary necessarily has a south A-edge adjacent 
to its left, and so ends up forming an A-labelled inner corner, as required for the labeling of the 
A-edges. Furthermore, it is clear that we obtain a partition after inserting the A-columns due to 
the interpretation of the lengths of the E-columns and the A-columns. 


It remains to show that the three conditions below are true for the multi-Catalan tableau if and only 
if they are true for the RAT it was mapped from. These are easily verified by carefully considering 
the rules for the fillings of both the RAT and the multi-Catalan tableaux, as seen in the example in 


Figure 15 


1. For any a, every tile above it in the same column is empty, there is no j3 to its right in the 
same row, and no a below it in the same column. 

2. Every /3 ends up in a DE or DA box, every tile to its left in the same row is empty, and there 
is no [3 to its right in the same row, and no a below it in the same column. 

3. Any empty box must have either an a below it in the same column or a /? to its right in the 
same row. 


We have thus a weight-preserving bijection from a RAT with the minimal tiling to a multi-Catalan 
tableau, and since all tilings of RAT are in bijection with each other, the Proposition follows. □ 


Finally, for the multi-Catalan tableaux, which translate to results for RAT with q = 0, there are 
some more refined enumerative results from |7j and also mm- 

Theorem 4.6. The weight generating function for the multi-Catalan tableaux of size n and whose 
type has r A’s is 




{anr-'-Y, 


2r p f2n — p\ a 


—^ 2n — p\n-\- r 
1=1 ^ 


v-p-l 


a 


-1 


/3-P-1 
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Theorem 4.7. The number of multi-Catalan tableaux of size n and whose type has r A’s is 


KrihhO) 


2(r + l) /2n + l 
n + r + 2\ n — r 


Theorem 4.8. Let n '■= r -\- k + 
has r A’s and k D’s is 


The number of multi-Catalan tableaux of size n and whose type 

r + l/n + l\/n + l\ 
n + l\ k y \ i )' 
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